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METHODS AND APPARATUS FOR THE SYSTEMIC CONTROL OF 
VENTILATORY SUPPORT IN THE PRESENCE OF 
RESPIRATORY INSUFFICIENCY 

5 This application claims the benefit of United States provisional application No, 

60/525,219, filed November 26, 2003. 

Field of the Invention 

This invention relates to methods and apparatus for the provision of ventilatory 
assistance synchronized with a patient's respiratory cycle and controlled to provide a 
10 minimum sufficient level of support in the presence of respiratory insufficiency. 

Background of the Invention 

Respiratory insufficiency is an abnormality of the respiratory system with an 
adverse effect on gas exchange, resulting in a lowering of the arterial partial pressure of 
oxygen, an increase in the arterial partial pressure of carbon dioxide, or both. It may be 
15 due to a decrease in ventilation in the presence of normal lungs, a lung abnormality, or a 
combination of these. 

Respiratory insufficiency generally results in worse gas exchange during sleep, 
for a variety of reasons including decreased respiratory drive and postural factors. The 
decrease in voluntary muscle function during REM sleep may cause a marked 
20 worsening of respiratory function during this stage of sleep, depending on the 
importance of the accessory muscles of respiration in the particular patient. 

Sleep may be associated with varying degrees of upper airway obstruction, 
referred to as obstructive sleep apnoea (OSA). 

Sleep-disordered breathing (SDB) generally refers to types of breathing 
25 disruption that occur during sleep. The most common form of sleep-disordered 

breathing is obstructive sleep apnea (OSA). Loud, intermittent snoring, apneas, and 
hypopneas characterize OSA. 
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Respiratory insufficiency, and possibly OSA, occur in conjunction with conditions 
experienced by patients with chest wall, neuromuscular, amyotrophic lateral sclerosis 
(ALS), or lung disease, such as chronic obstructive pulmonary disease (COPD). 
Because the symptoms of sleep apnea present themselves as a result of a precursor, 
5 SDB has become the general term used to describe any disease state that manifests 
apneas and/or hypopneas during sleep. Apneas and hypopneas interfere with gas 
exchange, fragment sleep, and frequently cause oxygen desaturations. In severe cases, 
patients may experience these oxygen desaturations and arousals from sleep hundreds 
of times each night. 

10 The most common treatment of OSA is to administer continuous positive airway 

pressure (CPAP). CPAP was invented by Sullivan and taught in U.S. Patent 4,944,310. 
Briefly stated, CPAP treatment acts as a pneumatic splint of the airway by the provision 
of a positive pressure, usually in the range 4-20 cm H 2 0. The air is supplied to the 
airway by a motor driven blower whose outlet is coupled by an air delivery hose to a 

15 nose (or nose and/or mouth) mask sealed with the patient's face. An exhaust port is 
provided in the delivery tube proximate to the mask. 

Ventilatory assistance may be provided by bMevel ventilators, proportional assist 
ventilators and servo-controlled ventilators. Each type of ventilator utilizes different 
methods for assisting with patient respiration and achieves different goals. 

20 Such ventilator devices provide appropriate responses to the changing conditions 

of the patient. For example, ventilatory devices determine when to trigger and cycle 
varying pressure levels associated with inspiratory and expiratory support so that the 
device will synchronize with the respiratory cycle of the patient. Triggering is the event 
associated with the initiation of the pressure levels intended for the patient's inspiration. 

25 Cycling is the event associated with switching to the pressure levels intended for the 
patient's expiration. Also, the devices may provide some method for increasing or 
decreasing ventilation during periods of hypoventilation or hyperventilation respectively. 
Maximizing machine performance in either or both of these areas generally results in 
greater patient comfort and better treatment of respiratory insufficiency 

30 Simple bi-level ventilators provide a higher pressure during the inspiratory portion 

of the patient's breathing cycle, a so-called IPAP, and a lower pressure during the 
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expiratory portion of the breathing cycle, a so-called EPAP. Traditionally, the switching 
may be accomplished by monitoring the respiratory flow or pressure and defining a 
threshold level. When the measured value exceeds the threshold, the device will trigger 
the IPAP pressure. When the measured value falls below the threshold, the device will 
5 cycle to the EPAP pressure. Other alternatives to such switching involve recorded 

respiration rates and the monitoring of elapsed time from the start of either inspiration or 
expiration. The machine may switch to the following portion of the respiratory cycle, 
either inspiration or expiration, after reaching the expected time for the previous part of 
the respiratory cycle. 

10 A different servo-ventilator device developed by ResMed Ltd. accomplishes 

synchronization by delivering smooth cyclical pressure changes based on a calculated 
instantaneous phase. Embodiments of the apparatus are the subjects of commonly 
owned U.S. Patent Application Serial No. 09/661,998 and Patent Nos. 6,532,957 and 
6,532,959, the disclosures of which are incorporated herein by reference. Generally, the 

15 apparatus provides an instantaneous mask pressure P(t) based upon a fraction of the 
patient's airway resistance R (this fraction ranging from zero to a substantial value but 
less than 1), respiratory airflow f(t), an amplitude A, and an estimation of the patient's 
instantaneous respiratory phase O as applied to a pressure waveform template tt(0) as 
follows: 

20 P(t) = P 0 + R f(t) + A tt(4>) for all f(t) (inspiration and expiration) 

where: 

A = Gj(V(t)-V TG T)dt 

Amin < A < Amax, and 

P 0 is an initial pressure. 

25 In this type of ventilator, V(t) can be, for example, one half the absolute value of the 
respiratory airflow f(t). The ventilation target V T gt may be a percentage of a measured 
volume of airflow, e.g., 95% of the average minute volume or a preset prescribed minute 
volume. G is the gain of the integral servo-controller, values in the range 0.1 - 0.3 
cmH 2 0 per L/min error in ventilation per second being suitable. Amin and Amax are limits 
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set on the degree of support A for comfort and safety. The limits 0.0 and 20.0 cmH 2 0 
respectively are generally suitable. 

In detecting the patient's respiratory phase, the apparatus uses a respiratory 
airflow signal and its derivative as input data for a set of fuzzy logic rules that are 
5 associated with particular phases of respiration. Using the results of the evaluations of 
the rules, a single continuous phase variable is derived and used as the instantaneous 
respiratory phase. This phase value as applied to a pressure waveform template then 
proportionally varies the delivered pressure in a manner that generates a realistic and 
comfortable respiratory cycle. Simultaneously, the calculation of A based on the target 
10 ventilation V T gt guarantees a desired level of ventilation. 

With regard to the issue of synchronization, while determining continuous 
instantaneous phase has its benefits in detecting the patient's respiratory cycle, it may 
be desirable to use a method for determining the pressure to be delivered to a patient 
which does not use a simple function n(^) to determine the number by which A in the 

15 above equation is multiplied. It may also be beneficial to utilize alternative methods for 
triggering the initiation of inspiratory pressure. 

With regard to the control of a pressure response to changing patient respiratory 
needs, in the above equation the rate of change of pressure support is simply 
proportional to the difference between the target ventilation and the actual ventilation. It 

20 may be desirable for small errors to result in a gentler response, and in particular for 

marked hypoventilation to result in a more brisk response, since it is much more likely to 
result in significant hypoxia. By contrast, when measured ventilation is well above 
target, it may not be necessary or desirable to briskly decrease the ventilatory support 
level, particularly considering that when there has been a sudden change in leak in a 

25 noninvasive ventilatory system, the measured ventilation is almost always greater than 
the actual ventilation until the leak estimation system has substantially compensated for 
the change in leak. 

A problem with ventilators that target total ventilation (traditionally called "minute 
ventilation") is that total ventilation may not provide an accurate measure of the extent to 
30 which the actual needs of the patient are met. Such a measure does not account for the 
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fact that patients have anatomical and physiological deadspace, and that the deadspace 
varies between patients. In particular, with a particular total ventilation, a high respiratory 
rate and a low tidal volume will provide lower alveolar ventilation than a low respiratory 
rate and a high tidal volume. A device which servocontrols minute ventilation is thus at 
5 risk of providing inadequate alveolar ventilation at high respiratory rates. 

Brief Description of the Invention 

It is an objective of the invention to provide a ventilator synchronized with a 
patient's respiratory cycle to provide a higher level of pressure support during inspiration 
and a lower level of pressure support during expiration. 

10 It is a further objective of the invention to provide such a synchronized ventilator 

that is flow triggered but phase cycled. 

It is a still further objective of the invention to provide a ventilator that maintains 
at least a target ventilation and responds to the error in ventilation in a way which is not 
linear with respect to the error and not symmetrical about zero error. 

15 It is another objective of the invention to provide a ventilator which uses two 

measures of ventilation which respond at different speeds to allow both smooth variation 
in pressure support and gain being a function of error, while still providing rapid 
response to marked hypoventilation. 

It is another objective of the invention to provide a ventilator which servocontrols 
20 alveolar ventilation. 

Additional objectives will be apparent to those skilled in the art upon 
consideration of the following description of the invention. 

One aspect of the invention involves patient-ventilator synchronization. In 
accordance with the invention, flow is measured. With the measure of flow, an 
25 instantaneous breath phase variable is determined using a fuzzy logic method. The 

ventilator then will trigger an inspiratory pressure for ventilatory support intended for the 
patient's inspiration when the measured flow exceeds a threshold. The ventilator 
subsequently cycles to deliver an expiratory pressure when the instantaneous phase 
passes a certain threshold value, such as 50%. 
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Another aspect of the invention relates to the delivery of servo-controlled 
ventilatory support to maintain a target ventilation. Adjustment of the delivered 
ventilatory support is based on the error between a target and actual measure, and a 
gain which varies with the distance of the actual measure from the target. Alternatively, 
5 it may be a function of relative ventilation error (which is defined to be ventilation error 
divided by target ventilation). This varying gain permits subdued ventilatory support 
responses for mild changes in patient ventilation while permitting more vigorous 
responses to more severe changes. In a simple form, a step function provides a lower 
gain if the distance from the target is less than a threshold and higher gain otherwise. 

10 The gain may be altered symmetrically around the target or alternatively asymmetrically. 
In an asymmetrical embodiment, the gain may change differently depending on whether 
the patient's ventilation is above or below the desired or target ventilation. These 
tempered or graduated gain adjustments permit the device to make a relatively greater 
increase (i.e., a more aggressive response) in pressure support when measured 

15 ventilation is further from the target, i.e., there is a large target error, than when 

measured ventilation is closer to the target, i.e., there is a small target error. It also 
permits making less aggressive changes to the gain for hyperventilation and more 
aggressive changes for hypoventilation. 

In another aspect of the invention, control of ventilation is more accurately 
20 achieved by controlling a measure of alveolar ventilation. Since the airways contain 
dead space, only a portion of the air inhaled by a patient can participate in gas 
exchange. In accordance with the invention, the total ventilation is corrected by 
subtracting the effect of anatomical and/or physiological dead space as an estimate of a 
target alveolar ventilation. This alveolar ventilation is utilized as a target and compared 
25 with a similarly calculated actual measure made during the delivery of ventilatory 
support. The target error of these two terms is then applied to adjust the pressure 
support provided by the ventilator. 

According to the invention, an error signal is computed that is the difference 
between a function of respiratory airflow over a period of time and a target value. Using 
30 a servo loop, breathing gas is delivered to the patient at a pressure that is a function of 
the error signal, the phase of the current breathing cycle, and a loop gain that varies 
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depending on the magnitude of the error signal. The loop gain increases with the 
magnitude of the error signal, and the gain is greater for error signals below a ventilation 
target value than for error signals above the ventilation target value. The target value is 
an alveolar ventilation that takes into account the patient's physiologic dead space. 

5 Brief Description of the Drawings 

FIG. 1 depicts the structure of the preferred ventilator apparatus for implementing 
the methods of the current invention; 

FIG. 2 is a flow chart of the general steps in synchronizing and setting the 
pressure; 

10 FIG. 3 is a flow chart showing steps to achieve flow triggering and phase cycled 

machine synchronization; 

FIG. 4 is a graph of Measured Ventilation V(t) versus Gain G illustrating a 
symmetrical function for adjusting the servo-ventilator gain that is a function of the 
distance that the measured ventilation deviates from the target ventilation; 

15 FIG. 4A is a graph of relative ventilation error versus gain illustrating an 

asymmetrical function for adjusting the servo-ventilator gain that is a function of the 
distance that the measured ventilation deviates from the target ventilation; 

FIG. 5 is a flow chart illustrating steps in the delivery of pressure in accordance 
with a target alveolar ventilation; and 

20 FIG. 6 is a graph of a function for determining dead space as a function of tidal 

volume. 

Detailed Description of the Invention 

In reference to FIG. 1, the ventilator relies on a servo-controlled loop for 
governing the operation of blower 2. The device includes a flow sensor 4f, pressure 
25 sensor 4p, a mask 6, and an air delivery conduit 8 for connection between the blower 2 
and the mask 6. Exhaust gas is vented via exhaust 13. Mask flow may be measured 
using a pneumotachograph and differential pressure transducer to derive a flow signal 
F(t). Alternatively, the pneumotachograph may be replaced by a bundle of small tubes 
aligned in parallel with the flow, from the blower with the pressure difference measured 



7 



WO 2005/051469 



PCT/AU2004/001651 



by the differential pressure transducer across the bundle. Mask pressure is preferably 
measured at a pressure tap using a pressure transducer to derive a pressure signal 
Pmask(t). The pressure sensor 4f and flow sensor 4p have been shown only symbolically 
in FIG. 1 since those skilled in the art will understand how to measure flow and pressure. 
5 Flow F(t) and pressure P maS k(t) signals are sent to a controller or microprocessor 15 to 
derive a pressure request signal P R e q uest(t). The controller or processor includes 
integrated chips, a memory and/or other instruction or data storage medium to 
implement the control methodology described in more detail herein. Programmed 
instructions with the control methodology are either coded on integrated chips in the 
1 0 memory of the device or loaded as software. 

A. Pressure Delivery Synchronization - Triggering and Cycling 

One aspect of the invention involves patient ventilator synchronization. In the 
preferred device, synchronization is performed by a unique combination of two different 
respiratory cycle detection techniques that collectively result in superior synchronization. 
15 It has been found that a ventilator controlled by application of flow-based triggering in 
conjunction with phase-based cycling as described herein has a synergistic effect in 
terms of improved machine synchronization. 

In general, the ventilator delivers a higher pressure upon the detection of patient 
inspiration and a lower pressure during expiration. One simple form of the delivered 
20 pressure as cycled by the controller is as follows: 



P(t) 



= P 0 + A 



if inspiration is true 



P(t) 



if inspiration is not true (i.e. expiration) 



where: 



A = 



Gj(V(t)-V TGT )dt 



25 



0 < 4^ < A < (chosen for safety and comfort) 



P 0 is an initial pressure, and 



G is a gain. 
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While the above pressure formula will generally result in a higher pressure and a 
lower pressure during the Inspiratory and expiratory portions of a respiratory cycle of a 
patient in a bi-level sense, those skilled in the art will recognize that the pressure 
changes may be ramped to provide more comfortable transitions from the higher to 
5 lower pressure or lower to higher pressure such that the ramping simulates a more 
natural respiratory airflow cycle. 

The steps in the methodology for triggering to the higher pressure and cycling to 
the lower pressure are listed in the flow chart of FIG. 2. In step 20, the current state of 
the machine is initialized to inspiration. In step 22, an instantaneous airflow is 

10 measured. In step 24, the patient's instantaneous phase is calculated. In step 26, the 
setting of the machine's state (inspiration or not inspiration) is based upon an evaluation 
of the airflow and instantaneous phase to determine the respiratory state of the patient. 
In step 28, the pressure is adjusted based upon the detected respiratory state. The 
steps repeat at a very high frequency relative to a normal respiratory cycle with the 

15 exception of the initialization step 20. 

With regard to steps 24 and 26, those skilled in the art will recognize an 
appropriate threshold and instantaneous flow measure suitable for detecting the onset of 
inspiration, which will also be adjusted to alleviate the effect of leak. As it relates to the 
calculation of instantaneous phase, the invention contemplates a continuous phase 
20 variable as disclosed in U.S. Patent Application Serial No. 09/661,998 or Patent No. 
6,532,957. The instantaneous phase variable is calculated in a fuzzy logic process. 

More specifically, the device utilizes a respiratory airflow and preferably its rate of 
change. The degree of membership of the respiratory airflow signal in each of the fuzzy 
sets zero, positive, large positive, negative, and large negative is calculated using 

25 suitable membership functions. Similarly, the degree of membership of the derivative of 
the respiratory airflow signal in each of the fuzzy sets steady, increasing, increasing fast, 
decreasing and decreasing fast is calculated using suitable membership functions. The 
degrees of membership in these sets are used in a group of fuzzy logic inference rules. 
Each fuzzy inference rule is associated with a particular phase of respiration. Although 

30 many variations are possible, in the preferred embodiment the inference rules relating to 
respiratory airflow are as follows: 
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1. If the airflow is zero and increasing fast, then the phase is 0 revolutions. 

2. If the airflow is large positive and steady, then the phase is 0.25 revolutions. 

3. If the airflow is zero and falling fast, then the phase is 0.5 revolutions. 

4. If the airflow is large negative and steady, then the phase is 0.75 revolutions. 

5. If the airflow is zero and steady and the 5-second low-pass filtered absolute 
value of the respiratory airflow is large, then the phase is 0.9 revolutions. 

6. If the airflow is positive and the phase is expiratory, then the phase is 0.1 
revolutions. 

7. If the airflow is negative and the phase is inspiratory, then the phase is 0.6 
revolutions. 

8. If the 5-second low-pass filtered absolute value of the respiratory airflow is 
small, then the phase in the respiratory cycle is increasing at a fixed rate equal to the 
patient's expected respiratory rate. 

9. If the 5-second low-pass filtered absolute value of the respiratory airflow is 
large, then the phase in the respiratory cycle is increasing at a steady rate equal to the 
existing rate of change of phase, low-pass filtered with a time constant of 20 seconds. 

The general method for developing and using these fuzzy inference rules is 
described in U.S. Patent Application Serial No. 08/935,785 or Patent No. 6,532,957. 
Generally, various features, such as the point of start of inspiration, are identified on a 
graph of airflow versus phase, and for each phase corresponding fuzzy rules are 
developed. For example, a suitable rule for the point "start of inspiration" could be "effort 
signal is small and the second derivative of the effort signal with respect to time is large 
positive." Membership functions, would cause that rule to be maximally activated at or 
near the start of inspiration. Preferably, the exact phase at the moment of maximal 
activation should be determined empirically. In the current example, the maximum 
activation will be at a phase shortly after the actual moment of start of inspiration, say 
0.05 revolutions, and this is the best phase to associate with the rule. The more features 
that are identified and assigned a rule and a phase, the smoother will be the resultant 
determination of instantaneous phase. Additional rules may be supplied so that the 
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phase will continue to rotate at an expected respiratory rate in the presence of an 
apnea or leak. Sample rules 8 and 9 serve these purposes respectively. 

Each of the rules in the combined set of fuzzy inference rules will be evaluated to 
determine a degree of activation G(n) by using a standard fuzzy inference method. For 
5 example, with respect to rule 2, using one such method assuming a unit weighting of 
rules, if (a) the degree of truth for the membership function "the airflow is large positive" 
evaluates to 0.6 and (b) the degree of truth for the membership function "the airflow 
signal is steady" evaluates to 0.4, and a fuzzy logic "AND" operator is applied, then the 
degree of activation for G(2) would be 0.4. 

10 Additionally, each of the fuzzy inference rules associates a particular rule with a 

particular phase O(n). For example, as shown above, rule 2 is associated with <£(2) = 
0.25 revolutions. Then, using the degree of activation G(n) for each <I>(n), a single value 
representing the instantaneous respiratory phase <b is calculated in a defuzzification step 
using the formula: 

15 <D = arctan ( 2 [G(n) sin O(n)] , Z [G(n) cos <£(n)] ). 

where the arctan function is a four-quadrant function, giving a phase of 0 in the unlikely 
event that both its arguments are zero. In the preferred embodiment, the instantaneous 
respiratory phase O is a continuous variable between 0 and 1 revolution, which 
proportionally represents a position within a singe respiratory cycle (starting at the 
20 beginning of inspiration (0) and ending at the conclusion of expiration (1)). Of course, 
those skilled in the art will recognize that the instantaneous respiratory phase can be 
represented by a different numerical range. 

As previously described, the delivered pressure is flow-triggered and phase 
cycled. This change in machine state is accomplished in step 26. One algorithm for 
25 achieving the invention is as follows: 

If (inspiration is true) and (instantaneous phase is. greater than an 
instantaneous phase threshold), then: 

set Inspiration to false (i.e., expiration detected) 
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Otherwise, if (inspiration is false) and (instantaneous airflow is greater than a 
flow threshold), then: 

set Inspiration to true (i.e., inspiration detected). 

The steps in such an algorithm are illustrated in the flow chart of FIG. 3. In step 
5 30 the apparatus examines the current respiratory state. If the current state is 

inspiration, the instantaneous phase is compared to a phase threshold in step 32. The 
phase threshold in the preferred embodiment is 0.5 revolutions but would be different 
depending on the numerical range of the instantaneous phase as it has been associated 
with one compete respiratory cycle. A central number in the range that would represent 
10 the switching point between inspiration and expiration would be appropriate. In step 32, 
if the instantaneous phase is greater than the threshold, step 34 is executed, otherwise 
no change to the machine's state will occur. In step 34, the machine's state will be set to 
expiration. 

If the state is not inspiration in step 30, step 36 compares an instantaneous 
15 measured flow with an airflow threshold. If the flow is greater than the threshold, a 

change of machine state will occur in step 38, otherwise no change will occur. In step 
38 the machine's state will be switched to inspiration. 

With this process, when the measure of the instantaneous flow exceeds the flow 
threshold while the machine is delivering an expiratory level of pressure, the machine 

20 will trigger the inspiratory pressure because the apparatus has effectively detected the 
onset of inspiration. Moreover, when the calculated instantaneous phase exceeds the 
phase threshold while the machine is delivering an inspiratory level of pressure, the 
machine will cycle to the expiratory pressure because the apparatus will have effectively 
detected the onset of expiration or otherwise detected a need to cycle to the onset of 

25 expiration. 

B. Controller Gain Varying with Ventilation Error 

In order to insure a proper ventilation of a patient in the presence of varying 
respiratory effort and patient condition, the pressure delivered P(t) during inspiration is 
set as a function of a target ventilation V^et and actual ventilation Vactuai, by adjustment 
30 of amplitude A. With a feedback measure of ventilation V actU ai, the support is servo- 
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controlled to satisfy the target ventilation V^et- The calculation of the delivered 
pressure support amplitude A is as follows: 

A-jG(B(t)).B(t)dt 

where E is a measure of ventilation error, such as absolute ventilation error: 

or relative ventilation error, defined by 

^target 

It is trivially the case that for some function H 
G(E{t)).E{t)=H(E{t)) 

10 The equation above is expressed in terms of G rather than H to indicate that in 

dA 

this integral servocontroller, — generally increases with E , but that this increase is 

at 

modified by the factor G{E) . 

Typically, when E is absolute ventilation error, a gain G of the integral servo- 
controller of a value in the range 0.1 to 0.3 cmH 2 0 per L/min error in ventilation per 

15 second is suitable. However, in servo-controlling the level of ventilation around the 
target, rapidly and completely meeting the target is not always desirable or necessary. 
For example, in some short instances of moderate hypoventilation experienced by the 
patient (i.e., V ac tuai falls below Vtarget), it may not be desirable to make a complete and 
immediate compensation response to the ventilation shortfall. A more rapid or 

20 aggressive adjustment might arouse a sleeping patient or otherwise have an 

unsatisfactory impact on the synchronization of the machine with the patient. Such a 
change may also be less comfortable for the patient. Accordingly, in the preferred 
embodiment of the invention, the gain G varies depending on the amount by which the 
measure of ventilation deviates from the target. This permits a slower response to the 

25 error when a need for a change in delivered ventilation is detected for minor deviations 
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from the target while also providing a more rapid response to the error for greater 
deviations from the target. In other words, the gain allows a patient to vary his 
ventilation around a target ventilation without an aggressive response if he does not 
stray too far. But if the ventilation does deviate from the target more significantly, a more 
5 aggressive response may be taken. Similarly, the gains may be distinct depending on 
whether the patient's ventilation is above or below the target so that the gain will allow 
the ventilation support to withdraw more gradually when the patient is over-ventilated 
since it might not be desirable to rapidly withdraw support that has recently increased as 
it might cause the patient to oscillate back into being under-ventilated, and the ventilation 
10 may appear to be above target because the leak has changed rapidly and the leak 
compensation system has not substantially compensated for this change. 

Various methods may be utilized to vary the gain as a function of the distance of 
the ventilation measure from the ventilation target. In one embodiment as illustrated by 
symmetrical step function 40 in the graph of FIG. 4, the gain G is a function of a 
15 threshold (T) 42. A higher gain per error is implemented for distances greater than the 
threshold and a lower gain per error for distances below the threshold. For example, the 
following algorithm can achieve this embodiment: 

Ifd<Tthen 

G = Lower Gain per error 

20 Otherwise 

G = Higher Gain per error 

Where: 

d is the distance from the target |V ac tuai - Vtargetl and 

T is a threshold distance. 

25 For example, the typical gain servo of 0.3 cmH 2 0 per L/min per error may be 

utilized as the higher gain per error and some fraction of higher gain per error would be 
appropriate for the lower gain per error (e.g., 0.2 cm H 2 0). In one embodiment, the 
threshold is a 1.0 liter per minute error where the ventilation measure and target 
ventilation are minute volumes. 
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The function 44 is another example in which the gain is varied as a function of 
the distance of the ventilation measure from the ventilation target. This function is also 
symmetrical around the target but does not exhibit sharp changes in gain. Moreover, 
unlike the symmetrical response shown in the previous embodiment, the gain 
5 adjustment may also be implemented asymmetrically to permit different gain responses 
depending on whether the patient's measured ventilation is above or below the 
ventilation target. These principles are illustrated by the preferred function shown in the 
graph of Fig. 4A. Those skilled in the art will recognize how to implement the graphed 
function of relative ventilation error with a ventilator device as illustrated in Fig. 4A. 

10 As shown in Fig. 4A f as measured ventilation moves away from the target, the 

gain may be graduated to increase the ventilation response of the ventilator. However, 
different gains are applied below and above the target. For respiratory airflows that are 
above the target value the gains are less than the gains for corresponding respiratory 
airflows that are below the target value. Thus, for a patient hyperventilating with the 

15 relative ventilation error being in excess of 0.5, the gain will not exceed some 

intermediate gain, e.g., 0.1 cm H 2 0/(L/min)/sec, since a more aggressive response is 
not critical. Similarly, since a leak in the system will change the patient's measure of 
ventilation to make it appear that the patient is hyperventilating, an aggressive response 
may be incorrect. Conversely, when the patient is hypoventilating to a significant degree 

20 (e.g., a relative ventilation error of -0.3), a greater gain is utilized, e.g., 0.3 cm 

H 2 0/(Lymin)/sec., since a more aggressive response is preferred. In the preferred 
embodiment, the measure of ventilation is a minute volume or minute ventilation; thus, 
the gain adjustment is a function of the relative minute ventilation error. 

When a gain which varies with error is used, the measure of ventilation used 
25 should preferably not change markedly during a breath. Using the measure of 

ventilation given in the above equation, namely, half the absolute value of respiratory 
airflow, the change in pressure support during the breath is the integral of the product of 
the gains and instantaneous errors for a range of ventilations, from a ventilation of zero 
at the transitions between inspiration and expiration to values of ventilation well above 
30 target ventilation at the times of peak inspiratory and expiratory flow (considering a 
typical breath when ventilation is at target). This works satisfactorily with a constant 
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gain, since gain can then be taken outside the integral, as in the equation above. But 
with a gain varying with error, such a measure of ventilation which changes substantially 
during the breath would not allow the desired gain for the ventilation typical of that breath 
to be delivered; instead the gains for a wide range of ventilations would be traversed, 
5 and in particular that for maximum ventilation error would be applied every time the 
respiratory flow crossed zero. Hence it is necessary to low-pass filter half the absolute 
value of respiratory flow in order to calculate ventilation and ventilation error before 
determining the gain. A disadvantage of this low-pass filtering is a delay in the response 
to changes in ventilation. 

1 0 One way of using a well-filtered ventilation measure for much of the time but 

being able to respond to severe hypoventilation quickly is to use two measures of 
ventilation with different speeds of response. The more rapidly responding measure of 
ventilation ("fast ventilation"), which fluctuates somewhat during a breath, has no effect 
on pressure support unless it is moderately below target, from which point as fast . 

15 ventilation decreases further the more slowly responding measure ("slow ventilation") 
becomes progressively less important and fast ventilation becomes progressively more 
important in the calculation of pressure support, until a point is reached below which fast 
ventilation is the only determiner of pressure support. 

In particular, one implementation, using slow and fast relative ventilation errors 
20 E s and E f respectively, uses the following equation for the rate of change of pressure 

support amplitude A : 



dA 
dt 



G S E S , when E f < 0.3 
G f E f when E f > 0.5 

aG s E s +(\~a)G f E f when 0.3 < E f < 0.5, where 

0.5 

a = — 

0.2 



Suitable gain values are Gf =0.7 and G s =0.5 cm H 2 0/sec. (The pressure support 
amplitude A discussed above can be determined by integrating dA/dt.) 
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C. Alveolar Ventilation Target 

Another aspect of the invention involves the control of the ventilator to maintain a 
ventilation target that more accurately corresponds to actual gas exchange in the 
patient's respiratory system. This aspect involves servo-control of a measure of alveolar 
5 ventilation. While the device as heretofore discussed may utilize a simple minute 
volume as the target ventilation, in this aspect of the invention the preferred target 
ventilation V TG t is an alveolar ventilation. 

FIG. 5 identifies steps in a methodology of servo-controlling the pressure to 
maintain alveolar ventilation. In an input step 50, a user is prompted to enter patient 
10 specific information necessary to assist with a calculation of an alveolar ventilation target 
for the patient. In step 52, an alveolar ventilation target is calculated. In step 54, a 
measure of alveolar ventilation is taken. In step 56, the measure and target are 
compared to derive an alveolar ventilation target error. In step 58, the pressure is 
adjusted based on the alveolar ventilation target error. 

15 In order for the device to control such a ventilation measure, an estimate of the 

patient's physiologic dead space V D (i.e., alveolar plus anatomical dead space) is 
required. In a simple embodiment, prior to the control of alveolar ventilation, as 
previously mentioned the user may be prompted to enter an estimate into the device in 
step 50. Thus, the user may enter alveolar volume V A or dead space volume V D 

20 determined from a prior analysis of the patient. V 0 may be determined by any known 
method for analyzing dead space such as by using Fowler's method or Bohr's method. 
The Fowler and Bohr methods for determining dead space volume are described by 
John B. West in Respiratory Physiology - the essentials 5^ ed., Baltimore MD, Williams 
& Wilkins, 1995. Alternatively, V D may be a fixed volume or it may be a default volume 

25 if no other information is entered or calculated, for example, 150 ml. 

Alternatively, the device may estimate anatomical dead space volume from a 
table based on user responses to questions regarding the size of the patient. For 
example, prior to use, the device may prompt a user for height, age and sex. With this 
information, the device may access a table or tables with such information associated 
30 with ideal weights for such a user. With the determined ideal weight of such a patient, 
anatomical dead space may be estimated by allowing 2 ml of anatomical dead space 
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volume for each 1 kg of idea! weight As those skilled in the art will recognize, the table 
may contain the anatomical dead space based on age, height and sex rather than the 
weight based on the one-to-one approximation noted above. The ideal weights for 
patients based on age, height and sex are known in the art. In one embodiment, the 
device may prompt the user for the actual weight and anatomical dead space volume or 
it may be approximated therefrom based on the above formula. This weight may also be 
adjusted based on the patient's body mass, such as by entry of a body mass index (BMI) 
or a skin fold thickness. Those skilled in the art will recognize various ways for utilizing a 
table or other formula or calculation to approximate anatomical dead space volume from 
the size of the patient. In one embodiment, the above weight-related estimate of 
anatomical dead space may be taken as V D or the estimate for both the alveolar and 
anatomic dead space. Alternatively, a separate estimate for alveolar dead space may 
be made and combined with the anatomical dead space. Thus, additional prompting of 
the device may also provide an estimate of alveolar dead space Vald that will be 
combined with the anatomical dead space Vand to serve as the estimate of total dead 
space volume V D (i.e., the physiological dead space \\Z D = V A ld + V a nd])- 

In one embodiment of the invention, the dead space V D may be taken to 
be a function of the patient's expected tidal volume or a recently calculated tidal volume 
measured during the provision of treatment as follows: V D = F(V T ). In such an 
embodiment, the function may be a fixed percent in a range of about 0.2 to 0.35 or 
roughly about 30% or about 1/3 of the tidal volume for normal patients at rest and higher 
for patients with a lung disease such as Chronic Obstructive Pulmonary Disease at rest. 
Given that certain patients with lung disease will have a greater degree of dead space 
than normal patients, in one embodiment the device may prompt input from a user to 
determine if the patient has a lung disease or COPD (for example, Emphysema or 
Chronic Bronchitis). The fixed percent may be set to a value at the higher end of the 
range for such patients, otherwise it will be set to the lower end of the range if the 
response is in the negative. Alternatively, an affine function based on tidal volume may 
be utilized to derive V D and, optionally this affine function may also depend on the 
disease condition of the patient as prompted by the device. For example, as illustrated 
in Fig. 6„an affine function 60 may be utilized as a function of tidal volume to determine 
dead space. Optionally, different functions may be selected depending on whether the 
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patient suffers from COPD, or any such particular disease. Thus a specified volume 
offset 62 associated with the patient's disease or the alveolar dead space (V A ld) may be 
applied to the patient's function of tidal volume as illustrated in Fig. 6. In still another 
embodiment of the invention, the V D may be estimated with use of an electronically 
5 integrated monitor to measure arterial partial pressure carbon dioxide PC0 2 or 
alternatively expired C0 2 in the mask. The C0 2 data may be utilized to derive an 
alveolar ventilation target during a learning period in accordance with Bohr's equation 
and then used for servo-controlling support to equal the target alveolar ventilation. 

As previously noted, the preferred target Vtarget is an alveolar minute volume. 
10 During the calculation step 52 of Fig. 5, the device may be configured or programmed 
with the following formula to calculate alvreolar minute volume target. 

Alveolar minute volume (V a tveoiar-minute) = V A x R 

where 

V A is the alveolar volume (V A = V T - V D ) 

15 V T is the patient's expected tidal volume 

V D is the patient's dead space volume, and 

R is the patient's expected respiratory rate (breaths/min). 

Alternatively, the target may be set as a function of this calculated Va^ar-minute. For 
example, the target may be set to 95% of this calculated Vaiveoiar-minute- 

20 For purposes of the expected values utilized in this determination, in step 50 the 

device may also prompt for the entry of the patient's expected tidal volume V T and 
expected respiratory rate R. However, such patient specific information may be 
measured by the device during an awake learning period or calculated based upon the 
entry of some of the above-identified information. For example, the respiratory rate may 

25 be taken as a median respiratory rate over the last 75% of an awake learning period. 
Optionally, where the target calculation is based on measures taken from an awake 
learning period, such as the respiratory rate, the calculation of the target alveolar minute 
volume may be adjusted to account for respiratory differences associated with sleep as 
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described in commonly owned U.S. Patent No. 6,644,312, the disclosure of which is 
incorporated herein by reference. 

A mean alveolar ventilation during part of a learning period can be calculated as 
follows. For each breath, there is a tidal volume V T , and if there are N breaths over a 
5 time T and the alveolar tidal volume is given as a function of tidal volume by the 
function a(V T ) , then the mean alveolar ventilation is given by 

=£«(»-») 

But when a median alveolar ventilation is required, which is desirable in order to reject 
outliers and artifacts, a breathwise median is inappropriate since it gives the same 

10 weight to short breaths (typically lower in tidal volume) as to longer breaths, and so 
typically does not give the same result as mean ventilation in cases where it should, 
such as a distribution of alveolar ventilation with respect to time which is symmetrical 
with respect to the mean. This is undesirable. A timewise median is more desirable. 
Mathematically, a piecewise constant function of time is constructed by assigning to the 

15 time period of each breath the alveolar ventilation (alveolar volume divided by breath 
duration) during that breath, and calculating the median of this function. In practice, an 
approximation to this can be efficiently calculated by forming a histogram of alveolar 
ventilation. For each breath, the alveolar ventilation is calculated, and the corresponding 
entry in the histogram is incremented by an amount proportional to the duration of the 

20 breath. The median of the histogram is then calculated in the usual way. 

The measurement of the actual delivered alveolar ventilation V actua i experienced 
by the patient in step 54 is based on the same formula as the target calculation with the 
exception that the tidal volume V T and the respiratory rate are recent measures 
preferably determined continuously from the measure of respiratory airflow from the 
25 patient. For example, the respiratory rate is the recent average respiratory rate or the 
number of breaths taken by the patient over the previous one-minute period. The tidal 
volume is the volume from the previous breath or average volume from the breaths 
taken during the previous one-minute period. Tidal volume may also be determined by 
dividing a measure of minute ventilation by respiratory frequency (VmJfresp). Preferably, 
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minute ventilation can be derived by filtering half the absolute value of respiratory 
airflow with a fourth-order low pass Bessel filter with a time constant selected to achieve 
the minute average. The patient's dead space volume V D as entered may remain 
constant or vary as a function of tidal volume as described above. 

Alternatively, the actual delivered alveolar ventilation V^^iar (or the target based 
on expected values) may be calculated according to the formula: 

f 

J resp 



Valveolar - ~ ^~ — V D 



fre 



J 



where: 

V to tai is a ventilation derived from airflow such as half the absolute value of 
10 respiratory airflow, low pass filtered with a comer frequency chosen to give a response 
appropriate for the application (for control, significantly faster than the servocontrolling 
response, to avoid oscillations; for target determination, possibly quite slow), 

V D is the dead space volume determined from any one of the methods 
previously described, and . 

15 fresp is the respiratory frequency (breaths per minute) which may be 

calculated as the number of revolutions per minute of the instantaneous phase variable 
as described previously. 

The preferred method of calculating alveolar ventilation is as follows: 

Suppose that the alveolar volume delivered as a function of volume V inhaled 
20 above FRC (functional residual capacity, the lung volume at the end of expiration) is 
given by U(V) . Then, both during inspiration and expiration, the following value is 
input to a lowpass filter whose output is the measure of alveolar ventilation: 

iigcrfjlel* 

Vo . J 

w here Q is respiratory flow as a function of time, the prime symbol indicates the 
25 derivative, and t 0 is the time of the start of the most recent inspiration or expiration. 
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10 



Intuitively, the half the absolute value of respiratory flow is multiplied by the rate at 
which the alveolar volume is increasing, relative to volume. For the special case of 
excluding only the anatomical deadspace V DA from the calculated alveolar ventilation, 

this is equivalent to calculating the volume since the beginning of the most recent 
inspiration or expiration, and when this volume is less than the deadspace, making the 
input to the lowpass filter zero, and otherwise making the input equal to half the 
absolute value of respiratory flow. This is because here we have 



so that 



when 0<V <V DA 



V DA , otherwise 



whenO<V<V DA 



otherwise 



The advantages of using both inspiratory and expiratory flow in the calculation of 
ventilation, rather than either inspiratory or expiratory flow alone, include lower sensitivity 
to baseline error in flow, and a smoother filter output. 

Once the actual and target alveolar ventilations are determined, delivered 
15 pressure support amplitude A can be adjusted based upon the error, for example, as 
follows: 

A = -G J (Vaiveolar-minute(t) - Vtarget ) dt 

The delivered pressure is then set in accordance with the amplitude A as 
previously discussed. 

20 D. Decreasing Pressure Support Adjustment When Leak is Changing 

Rapidly, 

Rapid change in leak causes the respiratory flow baseline to be transiently 
incorrect, which almost always has the effect of making calculated ventilation higher than 
actual ventilation, with the result that pressure support is either inappropriately reduced, 
25 or increased more slowly than it should be. 
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Patent No. 6,532,957 mentioned above, by Berthon-Jones, describes one 



method of decreasing this effect, using a fuzzy quantity called "jamming". 

In the. present invention, to decrease this effect, a different fuzzy measure R , 
where 0 < R < 1 , of the extent to which leak is changing rapidly, is used to decrease the 
5 rate of pressure support adjustment, by multiplying the rate of adjustment of pressure 

support amplitude — , as calculated above, by . To calculate R , the estimated 
dt 

respiratory flow is processed with a lowpass filter having a peak hold tendency, 
producing an estimate L of residual leak (defined to be leak which is not compensated 
for) which typically responds more rapidly that the leak estimate does, and so warns that 
10 the leak estimate is incorrect, though by itself fluctuating too rapidly to be directly useful 
as a leak estimate for the purposes of calculation of respiratory flow. When the absolute 
value of L is small (for example, less than the target ventilation), R is 0; when the 
absolute value of L is large (for example, twice the target ventilation), R is 1; and in 
between, it is interpolated linearly. 

15 An implementation of the filter which produces L is as follows. The first stage is a 

4-pole Bessel lowpass filter with a corner frequency of 3.2/60 Hz. The second stage is 
an instant-rise, variable time constant first order lowpass fall filter, where "rise" here 
means an increase in the direction of the previous output away from zero; so, for 
example, if the previous value of the output was negative, a "rise" means becoming 

20 more negative. A "fall" means movement towards zero. The time constant is longer 
when the flow is further from zero. Let the filtered respiratory flow (output of the first 
stage) relative to the target ventilation (both in the same flow units) be given by 




♦respiratory, filtered 

v 

r target 



Then the fall rate constant (reciprocal of the fall time constant) in units of 1/sec is 



25 



given by 
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1/0.5 when Q R < 0.75 

1/5 when Q R >2 

g*-0-75 



a/5+(l-cr)/0.5 otherwise, where 



1.25 



The purpose of the second stage is to cause the output of the second stage to be 
held closer to those values of the output of the first stage which are further away from 
zero. (The output of the first stage typically has the appearance of a sinusoid plus a 
5 slowly changing "DC value.) 

This method of compensating for residual leak has the advantage over the 
"jamming" described in Berthon-Jones Patent No. 6,532,957 in that it is able to detect 
residual leak which is not large enough to produce apparent inspiratory or expiratory 
periods much longer than normal (to produce these long periods generally requires a 
10 baseline shift large enough to make estimated respiratory flow positive throughout or 
during most of expiration, for example, in the case where leak has suddenly increased). 

Although the invention has been described with reference to various 
embodiments as described in this entire disclosure, it is to be understood that these 
embodiments are merely illustrative of an application of the various principles of the 
15 invention. Numerous modifications, in addition to the illustrative embodiments of the 
invention discussed herein, may be made and other arrangements may be devised 
without departing from the spirit and scope of the invention. 
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